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CONVAIR, SCIENTIFIC  RESEARCH  LABORATORY 
^  RESEARCH  NOTE  1 

STAGNATION  POLNT  HEAT  TRANSFER  FOR  HYPERSONIC  FLOW 
Mary  F.  Romig 


SUMMARY 

A  design  method  is  given  for  obtaining  stagnation  point  heat  transfer 
on  a  spherical  nose  in  hypersonic  flow.  The  heat  transfer  equation  incor¬ 
porates  real  gas  effects  through  use  of  dissociated  air  properties  and 
a  fitted  curve  for  the  stagnation  point  velocity  gradient.  Comparison 
with  numerical  solutions  of  the  compressible  dissociated  stagnation  point 
heat  transfer  indicates  agreement  within  -  5%.  The  derived  equation  for 
heat  transfer  is  a  function  only  of  free  stream  Mach  number,  free  stream 
pressure  and  body  nose  radius. 


SYMBOLS 

C  ■  constant  in  heat  transfer  equation  (Eq.  (1)  ) 
h  *  enthalpy 

i  ■  dimensionless  enthalpy  *  h/h 

O 

i*  *  uimensionless  reference  enthalpy  (Eq.  (3)  ) 

k  *  ratio  of  free  stream  density  to  density  behind  the  normal  shock 
M  ■  Mach  number 
p  *  pressure 

Pr  -  Prandtl  number  *  viscosity  x  specific  heat/thermal  conductivity 
q  ■  compressible  stagnation  heat  transfer  rate 
Q  *  the  ratio  q/ vC 
«■  radius  of  the  nose 
R  -  gas  constant 

s  *  surface  distance  around  nose 

T  *  temperature 

u  *  velocity 

u'  -  velocity  gradient  «  du/ds 
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Greek 


y  -  ratio  of  specific  heats 

X  ■  factor  in  heat  transfer  equation  (Eq.  (13)  ) 

H  ■  viscosity 

p  ■  density 

Subscripts 
oo  »  free  stream 
o  »  stagnation 

w  ■  wal 1 

B  -  base  (at  T  -  400*11) 

D 

1  -  incompressible 

Superscript 

*  ■  quantity  evaluated  at  reference  enthalpy 

AHALYSIS 

The  equations  presented  in  this  note  are  the  result  of  a  semi- 
empirical  approach  to  the  problem  of  determining  stagnation  point  heat 
transfer.  This  approach  arose  from  the  fact  that  no  closed  exact 
solution  exists  for  the  differential  equations  determining  the  compres¬ 
sible  flow  heat  transfer  at  the  stagnation  point.  It  has  long  been 
the  practice  (see  e.g.,  Ref.  1)  of  aerodynamic  is ts  to  modify  the  well- 
verified  incompressible  flow  boundary  layer  heat  transfer  equatlonsso 
they  could  be  used  for  compressible  flow. 

Consider  the  Incompressible  stagnation  point  heat  transfer  equation, 

2 

given  by  Slbulkln 

q  »  C  nT^T  J  uj‘  (h.  -  h  )  (1) 

1  w 

where  C  ■  0.763  for  Pr  ■  1. 

An  incompressible  flow  heat  transfer  equation  can  be  used  to  obtain 
compressible  flow  heat  transfer  if  it  is  possible  to  account  for  the 
effects  of  compressibility  on  the  air  property  values  (pp)  and  the 
stagnation  point  velocity  gradient  (uj).  In  the  discussion  following, 


-3- 


these  two  problems  are  treated  separately  since  the  velocity  gradient  may 
be  determined  without  reference  to  the  boundary  layer  equations. 


Various  temperature  levels  have  been  suggested  at  which  to  evaluate 
the  properties  pp.  Among  them  are  the  wall  temperature ,  the  stream 
temperature,  or  an  ave'--;^  temperature.  Tht.  approximation  to  flat 

plate  compressible  heat  transfer  has  been  obtained  by  evaluating  pp  at 
a  "reference  temperature",  which  is  a  function  of  local  Mach  number  and 
wall-to-stream  temperature  ratio.  A  complete  discussion  of  this  technique, 
beyond  the  scope  of  this  paper,  may  be  found  in  Ref.  1. 

It  has  been  well  established^  that  compressibility  effects  for  flat 
plate  laminar  flow  can  be  introduced  through  use  of  such  a  reference 
temperature.  It  is  assumed  that  this  method  will  apply  to  stagnation 
point  flow  as  well.  If  it  is  further  assumed  that  real  gas  effects  are 
included  through  use  of  a  reference  enthalpy,  rather  than  a  reference 
temperature,  and  the  proper  air  properties,  then  Eq.  (1)  may  be  used 
to  determine  stagnation  point  heat  transfer  in  hypersonic  flow.  It  will 
be  seen  later  that  there  is  some  evidence  to  support  these  assumptions. 


In  order  to  simplify  the  analysis  several  further  assumptions  will 
be  made.  These  may  be  modified  without  invalidating  the  following 
development . 

(1)  Newtonian  flow  equations  are  used  for  surface  pressure  and 
stagnation  point  velocity  gradient 

(2)  1  and  h,^  h^,  i.e.,  the  wall  is  cool  compared  to  the  stream 

(3)  Pr  -  1 

If  assumptions  (2)  and  (3)  are  changed,  the  analysis  presented  here 
can  be  applied  to  supersonic  flow  as  low  as  M  *3.  The  curve-fitted 

GO 

equations  shown  later  for  u,  and  p|i  also  permit  this  extension. 
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The  following  assumptions  are  less  subject  to  modification. 

(4)  Free  stream  temperature  T  ■  400*R  ■  T 

(5)  The  body  shape  In  the  stagnation  region  Is  spherical. 

Although  assumption  (4)  greatly  simplifies  the  final  equations  it  was 
chosen  primarily  so  that  the  normal  shock  solution  of  Ref.  6  could  be 
utilized  for  the  velocity  gradient.  It  is  not  known  at  the  present  time 
how  large  an  effect  the  change  of  initial  temperature  would  have  on  the 
normal  shock  solution. 

On  the  basis  of  these  assumptions,  the  compressible  stagnation  point 
heat  transfer  becomes 


q 


0.763 


(2) 


for  b/personlc  flow,  The  dimensionless  reference  enthalpy,  defined  in 
Ref.  1,  reduces  to 


i*  -  0.5  (i  +  i  ) 
o  w 

at  the  stagnation  point.  From  assumption  (2) 


i*  »  0.5  i 


(3) 


where  i  *  h/h  . 

D 

The  stagnation  enthalpy,  determined  from  free  stream  conditions,  is 


h  7  *  1  , 

£  - 1  +  -v-  « 

h  L  oo 


For  h  *  h  (assumption  (4)  and  M  ■>>  1, 

co  B  co 


7  -  1 


M 


(4) 


2 
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k 


and,  from  Eq.  (3),  for  7  »1.4, 

00 

1*  -  0.1  M2  (5) 

00 

so  that  Eq.  (2)  becomes  (for  -  96.8  BTU/lb) 

J 

q  »  0.763  x  96.8  x  0.2  nJ  p^*1  J  up  M2  (6) 

.  *> 

Examination  of  Eqs .  (5)  and  (6)  shows  it  would  be  desirable  to  express 
p*p*  (or  1*)  and  u*  as  functions  of  given  (free  stream)  conditions.  It 
is  possible  to  do  this  through  use  of  fitted  carves  and  assumptions  (l)-(5). 

The  Air  Properties 

For  moderate  supersonic  Mach  numbers  the  values  of  p*p*  to  be  used 
in  Eq,  (6)  are  comparatively  well-defined.  As  the  air  temperature  Increases 
with  speed,  however,  the  various  constituents  become  dissociated  and  com¬ 
plex  chemical  reactions  occur.  In  this  analysis  the  air  will  be  assumed 
to  be  in  dissociation  equilibrium.  The  density  and  viscosity  selected  for 
use  in  this  analysis  are  associated  with  the  most  recent  complete  set  of 

air  properties^.  While  the  thermal  properties  (p,  T,  p,  h)  have  subse- 

4 

quently  been  revised  the  transport  properties  have  yet  to  be  computed.  It 
is  suggested  that  when  new  transport  properties  become  available  they  be 
Incorporated  into  the  analysis  in  a  manner  similar  to  the  following  procedure. 

Fig.  1  shows  the  variation  of  pp/pgM^  with  enthalpy.  The  pressure 
dependency  of  density  is  Incorporated  into  the  base  density  by  the 
expression 


n  _  P(local) 

PB  RT,, 

D 

where  T  ■  400*R.  The  variation  of  pp/p  u  with  pressure  due  to  dissociation 
is  slight  (within  -  67.  from  0,1- i0  atmospheres)  so  the  curve  for  p  ■  ] 
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atmosphere  Is  plotted  In  Fig.  1.  A  good  approximation  to  the  curve  of 
Fig.  1  Is  given  by 

..,-.28 
pp  “  pb^b  (i) 

or,  In  this  particular  case,  from  Eq.  (5), 

P*P*  -  1.9  Pg^H/'56  (7) 

Nov,  the  above  equation  Is  still  a  function  of  local  flow  properties 
through  the  density  p^.  The  local  pressure,  in  this  case  the  stagna¬ 
tion  pressure,  can  be  removed  through  use  of  the  Newtonian  flow 
approximation  at  hypersonic  Mach  numbers  (Kef.  7,  p.265), 

p-p  +  7  M2  p  £27  K2p 

O  00  00  00  00  000000 


so  that 


P*|i* 


‘■’Vi 
*T»  "• 


m1.44 


(8) 


The  Velocity  Gradient 

The  stagnation  point  velocity  gradient  for  hypersonic  flow  on  a 
sphere  can  be  given  by  the  constant-density  Newtonian  flow  approximation 
developed  by  Li  and  Geiger 

”  Us  \(  k(2  -  k)1  (9) 

*N 

where  k  is  the  ratio  of  free  stream  density  to  density  behind  the  normal 
shock.  Values  of  'T  k(2  -  k)1  are  glvei  in  Fig.  2  for  the  normal  shock 
solution  of  Ref.  5,  which  Is  based  on  the  recent  air  properties  of  Ref.  4. 
Fig.  ?  shows  that  u,'  can  be  approximated  by 
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0.8 


u  (M  ) 

00  09 


-.232  0.8 


*N 


{V  RT  '  M 

*  00  00  c 


.768 


(10) 


»  1  | 

neglecting  the  variation  in  \J  k(2  -  k)  with  free  stream  pressure.  The 
curve  given  by  Eq.  (10)  also  agrees  with  the  low  speed  test  data  of  Ref.  5. 


Another  method  of  obtaining  uj;  for  hypersonic  flow  is  to  consider  use 
of  an  effective  7  rather  than  use  of  real  gas  air  properties.  Lees^ 
suggests  the  formula,  which  at  high  Mach  numbers  reduces  to 


u 


i  ~ 

•  ~ 


u 

00 


(U) 


where  y  Is  suggested  to  vary  between  1.1  -  1.2  for  high  temperature  flow. 
This  equation  Is  also  shown  in  Fig.  2  for  y  «  1.1,  1.2  and  1.4.  It  can 
be  seen  that  using  a  constant  value  of  y  restricts  the  formula  for  uj 
to  a  given  Mach  number  region.  For  this  reason  Eq.  (10),  which  is 
applicable  over  a  wide  range  of  Mach  numbers,  Is  recommended  for  general 
use.  In  specific  Instances  where  the  Mach  number  range  permits  use  of 
a  constant  7,  Eq.  (11)  is  recommended  on  the  basis  of  its  s  jnpllclty. 


Stagnation  Point  Heat  Transfer 


The  UBe  of  Eqs.  (8)  and  (10)  now  allow  the  heat  transfer,  Eq.  (6), 
to  be  expressed  solely  as  a  function  of  free  stream  Mach  number,  pressure 

and  nose  radius.  After  the  constants  are  evaluated,  with  T  -  T_  *»  400°R, 

00  0 

then 


q 


0.0145  M3*1 

00 


(12) 


The  restrictions  on  Eq.  (12)  are  that  M  »  1,  the  wall  Is  much  cooler  than 

00 

the  stagnation  point  air  and  free  stream  temperature  T  ■  400*R.  These 

00 

conditions  could  correspond  to  hypersonic  flight  in  the  Isothermal 
altitudes . 


/V 
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COMP ARISON  WITH  OTHER  SOLUTIONS 


Unfortunately,  we  cannot  present  an  experimental  comparison  for  Eq.  (12) 
since  existing  test  data  in  this  flight  regime  are  classified.  In  most  cases 
of  numerical  solutions  for  the  compressible  stagnation  point  the  determina¬ 
tion  of  pp.  and  uj  are  left  in  the  equation  as  parameters  to  be  determined 
by  the  reader.  Thus,  the  only  verification  cf  Eq.  (12)  which  can  be  offered 
is  to  pheck  the  validity  of  the  use  of  a  reference  enthalpy,  i.e.,  to  deter¬ 
mine  how  well  the  method  approximates  the  conditions  of  compressibility  and 
dissociation. 


The  reference  enthalpy  method  can  be  tested  in  this  way  by  inserting 
the  curve-fitted  equation  for  p*p*  (Eq.  (7)  )  into  the  heat  transfer 

g 

equation  (Eq.  (2)  )  and  comparing  with  the  numerical  solution  of  Mark  . 

This  solution  of  the  compressible  stagnation  point  boundary  layer  equations 
Include  use  of  the  dissociated  air  properties  of  Ref.  3.  In  Mark's  equation* 


where  X  .  a  function  of  the  shear  stresB,  was  determined  by  numerical 
Integration  cf  the  boundary  layer  equations.  Eq.  (2)  can  be  written 


ht 


'  PB^B 


The  terms  in  brackets  should  compare  and  will  indicate  the  effectiveness  of 
the  reference  enthalpy  in  approximating  compressible  heat  transfer.  The 
ratio 

q(l«)  ».052(MJ--28PrM 

q(Mark)  * 

is  plotted  in  Fig.  3  for  the  range  of  flight  velocities  covered  by  Ref.  8. 
The  agreement  is  within  *.57.  for  the  three  wail  temperature  values  given 
by  Mark.  This  tends  to  substantiate  the  assumptions  made  previously  that 
the  reference  enthalpy  method  car.  be  used  at  the  stagnation  point,  and 
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furtherraore  that  it  accounts  for  real  gas  effects  on  the  compressible 
boundary  layer. 


Another  calculation  *«diich  can  be  made  points  out  the  effects  of 
compressibility  and  dissociation  on  the  stagnation  point  boundary  layer 
heat  transfer.  A  plot  of  such  a  calculation  also  offers  a  simple  graphical 
method  of  obtaining  heat  transfer  from  a  standard  solution.  The  function 
Q  ■  q/  f  u-'  can  be  determined  from  Eqs .  (2)  and  (7)  for  compressible  flow. 
The  incompressible  Qi  -  q^/  \j  uj ' is  given  in  Eq.  (1),  evaluated  at  a  given 
wall  temperature  T  .  The  ratio 


3- 


Jl.9(M  )“*28 
1  00 

*  V»'Vi ' 


(14) 


ie  plotted  in  Fig.  4  for  -  2Q00*R.  It  indicates  the  effect  of  compressi 
bility  and  dissociation  on  the  boundary  layer.  Since  is  comparatively 
easy  to  evaluate,  Fig.  4  presents  a  rapid  method  of  obtaining  Q,  which  with 
Eq.  (10),  gives  the  value  of  compressible  heat  transfer. 

It  can  be  seen  in  Fig.  4  that  evaluation  of  the  air  properties  at 
the  wall  conditions  could  overestimate  heat  transfer  by  as  much  as  30Z 
at  Mach  20.  Since  Eq.  (1)  is  sometimes  used  to  obtain  engineering  estimate 
the  possible  error  induced  should  be  kept  in  mind. 
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